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Abstract. Based on the refined dynamical model pro- 
posed by us earlier for beamed 7-ray burst ejecta, we 
carry out detailed numerical procedure to study those 7- 
ray bursts with rapidly fading afterglows (i.e., ~ t~ 2 ). It is 
found that optical afterglows from GRB 970228, 980326, 
980519, 990123, 990510 and 991208 can be satisfactorily 
fitted if the 7-ray burst ejecta are highly collimated, with 
a universal initial half opening angle 6q ~ 0.1. The obvi- 
ous light curve break observed in GRB 990123 is due to 
the relativistic-Newtonian transition of the beamed ejecta, 
and the rapidly fading optical afterglows come from syn- 
chrotron emissions during the mildly relativistic and non- 
relativistic phases. We strongly suggest that the rapid fad- 
ing of afterglows currently observed in some 7-ray bursts 
is evidence for beaming in these cases. 
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1. Introduction 

The cosmological origin of 7-ray bursts (GRBs) has 
been well established due to recent discovery of multi- 
wavelength afterglows (Costa et al. 1997; Metzger et al. 
1997; Galama et al. 1997; Wijers, Rees & Meszaros 1997; 
Piran 1999). However, we are still far from resolving the 
puzzle of GRBs, because their "inner engines" are well hid- 
den from direct afterglow observations. Some GRBs local- 
ized by BeppoSAX satellite have implied isotropic energy 
release of more than 10 ergs (Kulkarni et al. 1998, 1999; 
Andersen et al. 1999; Harrison et al. 1999), which forces 
many theorists to deduce that GRB radiation must be 
highly collimated. Obviously, whether GRBs are beamed 
or not has become one of the most important problems 
that need to be solved urgently. 



In the literature, it is generally believed that afterglows 
from jetted GRB remnant are characterized by an obvious 
break in the light curve during the relativistic phase, due 
to both the jet edge effect (Panaitescu & Meszaros 1999; 
Kulkarni et al. 1999; Meszaros & Rees 1999) and the lat- 
eral expansion effect (Rhoads 1997, 1999). The breaking 
point is determined by 7 ~ 1/0, where 7 is the Lorentz fac- 
tor of the jet and 9 is the half opening angle. Recently we 
have developed a refined dynamical model that can cor- 
rectly describe the overall evolution of an ultra-relativistic 
jet to non-relativisitic phase with the expanding veloc- 
ity as small as 10~ 3 c (Huang et al. 2000a). Surprisingly 
enough, our detailed numerical results (Huang, Dai & Lu 
2000b) show that the break theoretically predicted in light 
curve does not appear during the relativistic phase, i.e., the 
time determined by 7 ~ 1/9 is not a breaking point. How- 
ever, an obvious break does appear within the relativistic- 
Newtonian transition region, the degree of which is found 
to be parameter dependent (Huang, Dai & Lu 2000b). 
Generally speaking, the Newtonian phase of jet evolution 
is characterized by a rapid decay of optical afterglows, 
with the power-law timing index a ~ 1.8 — 2.1. 

In practical observations, the power-law decay indices 
of afterglows from GRB 980326, 980519 and 991208 are 
anomalously large, a ~ 2.0 (Groot et al. 1998; Owens et al. 
1998; Castro-Tirado et al. 1999b), and optical light curves 
of GRB 990123 and 990510 even show obvious steepening 
at observing time t > 1 — 2d (Kulkarni et al. 1999; 
Harrison et al. 1999; Castro-Tirado et al. 1999a). Recently 
GRB 970228 was also reported to have a large index of 
a ~ 1.73 (Galama et al. 1999b). These phenomena have 
been widely regarded as evidence for beaming (Sari, Piran 
& Halpern 1999; Castro-Tirado et al. 1999a). The purpose 
of this Letter is to study these cases numerically, based 
on our refined beaming model (Huang et al. 2000a). It 
is found that optical afterglows from these GRBs can be 
easily reproduced, thus a jet model is strongly favored. 
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Fig. 1. Observed optical (R band) afterglow from GRB 990123 
(Galama et al. 1999a; Kulkarni et al. 1999; Fruchter et al. 1999; 
GCN 354). The solid line is our best fit to it by employing a 
jet model. Inset shows the evolution of 7 in our model 



2. Model 

The importance of non-relativistic expansion phase to 
our understanding of GRB afterglows has been stressed 
very early by Huang et al. (1998a, b). A refined generic 
model that is appropriate for both ultra-relativistic and 
non-relativistic isotropic blastwaves has been proposed 
(Huang, Dai & Lu 1999a, b). Very recently, a similarly 
refined dynamical model for beamed GRB ejecta was also 
developed (Huang et al. 2000a). Our calculation here will 
be based on this model. For completeness, we describe 
the model briefly here. For details please see Huang et al. 
(2000a). 

Let R be the distance from the burster in the burster 
frame, M e j be the initial ejecta mass, n be the particle 
number density of the surrounding interstellar medium 
(ISM), and m be the swept- up ISM mass. The evolution 
of the beamed ejecta is described by (Huang et al. 2000a): 



= Pnh + vV - 1), 



= 2ttR 2 (1 — cos6)nm p , 

dR 



dt R 

d7 _ 7 2 - 1 

dm M j + em + 2(1 — e)-fm ' 



c g=^- 1 >fr- 1 > l + 7(7-l) c2 ' 



(1) 
(2) 

(3) 
(4) 
(5) 



where (3 = A/7 2 — 1/7, m p is the proton mass, c s is the 
co-moving sound speed, e is the radiative efficiency (here 
we consider only adiabatic jets, for which e = 0), 7 « 
(47 + l)/(3 7 ) is the adiabatic index (Dai, Huang & Lu 
1999). 
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Fig. 2. Observed R band light curve of GRB 990510 and our 
model fit to it. Data points were taken from GCN 310, 313, 
315, 321, 323, 325, 329 and 332 (also see Stanek et al. 1999; 
Harrison et al. 1999). Inset shows the evolution of 7 in our 
model 



A strong blastwave will be generated due to the in- 
teraction of the jet and the ISM. Synchrotron radiation 
from the shock accelerated ISM electrons gives birth to 
afterglows (Sari, Piran & Narayan 1998; Vietri 1997a, 
b). As usual we assume that the magnetic energy den- 
sity in the co-moving frame is a fraction £{| of the to- 
tal thermal energy density (B' 2 /8n — £|e'), and that 
electrons carry a fraction £ e of the proton energy. This 
means that the minimum Lorentz factor of the random 
motion of electrons in the co-moving frame is 7 c . m i n = 
£ e (7 — l)m p (p — 2)/[m e (p — 1)] + 1, where p is the index 
characterizing the power-law energy distribution of elec- 
trons, and m e is the electron mass. Our model also takes 
the electron cooling (Sari, Piran & Narayan 1998) and the 
equal arrival time surface effect (Panaitescu & Meszaros 
1998) into account. 

The most important advantage of this model is that 
it is appropriate for both adiabatic and radiative jets, no 
matter whether their expansion is highly relativistic or 
Newtonian. 

3. Numerical Results 

Based on the model described above, Huang, Dai & Lu 
(2000b) have found that: (i) The optical light curve does 
not break during the relativistic phase, i.e., the time deter- 
mined by 7 ~ 1/9 is not a breaking point, (ii) An obvious 
break does appear within the relativistic-Newtonian tran- 
sition region, but its existence depends on parameters such 
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Fig. 3. Observed R band light curve of GRB 980519 and our 
model fit to it. Data points were taken from GCN 82, 83, 87, 
88, 91, 148, 149 (also see Halpern et al. 1999). Inset shows the 
evolution of 7 in our model 
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Fig. 4. Observed R band light curve of GRB 991208 and our 
model fit to it. Data points were taken from GCN 452, 454, 
456, 458, 461 and IAU Circ. 7332. Inset shows the evolution of 
7 in our model 



as £ e , n i ^o- Increase of any of them to a large enough 
value will make the break disappear, (iii) Generally speak- 
ing, the Newtonian phase of the jet evolution is charac- 
terized by a rapid decay of optical afterglows (a ~ 1.8 - 
2.1). 

In this section we use the model to study those GRBs 
whose optical afterglows decayed rapidly. They include 
GRB 970228, 980326, 980519, 990123, 990510 and 991208. 
We employ a standard Friedmann cosmology with Hq = 
65 km s" 1 Mpc" 1 , Q = 0.2, A = throughout. 

3.1. GRB 990123 

GRB 990123 was determined to be at a redshift z = 
1.6004 ±0.0008 (Kulkarni et al. 1999), corresponding to a 
luminosity distance £>l ~ 12 Gpc. An enormous isotropic 
7-ray energy release of E 1 « 3.4 x 10 54 ergs was estimated 
(Kulkarni et al. 1999). R band observations made be- 
tween 0.16 and 2.75 days after the GRB are well described 
by a power-law with index a. — 1.09 ± 0.05 (Fruchter 
et al. 1999), and the light curve steepens greatly when 
t > 4 d. Fig. 1 illustrates our best model fit to the 
R band light curve. We have taken the following ini- 
tial values and parameters: initial energy per solid angle 
{E/n) = 3.4 x 10 54 ergs/47r, 7o = 300, n = 10 3 cm" 3 , 
Co = 0.2, £| = 10~ 6 , O = 0.13, p = 2.2 and D h = 12 
Gpc. In our calculation, the jet enters the sub-relativistic 
phase at time t ~ 10 5 5 s, leading to an obvious break in 
the light curve. We see that observed data points spanning 
from t = 0.16 d to t = 59.5 d can be well fitted. 

3.2. GRB 990510 

GRB 990510 lies at a redshift z > 1.619 ±0.002 (Vreeswijk 
et al. 1999), corresponding to £>l > 12 Gpc. Adopt- 
ing £>l = 12 Gpc implies an isotropic energy of E 1 w 



2.9 x 10 53 ergs (Harrison et al. 1999). Fig. 2 illustrates R 
band afterglows reported in the literature, where the solid 
line is our model fit, taking: (E/Q,) = 2.9 x 10 53 ergs/47r, 
70 = 300, n = 1 cm" 3 , £ e = 0.2, £| = 0.01, O = 0.075, 
p = 2.6 and Z?l = 12 Gpc. In our calculation, the ejecta 
enters the sub-relativistic phase at t ~ 10 6 s. But due 
to relatively large values of £ c and the theoretical 
light curve peaks very late (Huang, Dai & Lu 2000b), 
£ P eak ~ 10 4 s, so that observed data points before 10 5 
s can be well fitted. For the same reason (also see Huang, 
Dai & Lu 2000b), theoretical flux decays rapidly during 
the mildly relativistic phase (i.e., from 10 5 s to 10 6 s), 
which just meets the observational requirements. 

3.3. GRB 980519 

GRB 980519 had a rapid fading in optical as well as 
in X-ray band. Its optical afterglow is consistent with 
f -2.05±o.04 (Halpern et al. 1999). No redshift was deter- 
mined. Here we assume a typical value of z ~ 1.4 for it, 
corresponding to Z?l ~ 10 Gpc. Then the BATSE mea- 
sured 7-ray fluence of (2.54 ± 0.41) x 10~ 5 ergs cm -2 
implies an isotropic energy of E 1 ~ 10 53 ergs. Taking 
(E/n) = 1 x 10 53 ergs/47r, 7o = 300, n = 100 cm- 3 , 
£ c = 0.1, £l = 0.01, O = 0.1, p = 2.6 and D L = 10 Gpc, 
we give our model fit to the R band light curve in Fig. 3. 
In our calculation, due to the relatively large ISM density, 
the jet becomes sub-relativistic at t <~ 10 5 s, so the theo- 
retical afterglow decays rapidly after t ~ 10 4 s. It is clear 
that the rapid fading of GRB 980519 can be reasonably 
explained in our model. 

34. GRB 991208 

GRB 991208 was localized by the Interplanetary Network. 
Beginning on 1999 Dec. 10.27 UT, the observed opti- 
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cal afterglow is consistent with a rapid decay of t~ 2 - 5 . 
A redshift of z fa 0.71 was measured (Dodonov et al. 
1999), corresponding to L\ ~ 4 Gpc. The implied 7- 
ray energy is £L ~ 10 53 ergs. We give our model fit to 
the R band afterglows in Fig. 4, where we have taken: 
(E/n) = 1 x 10 53 ergs/47r, 7o = 300, n = 10 cm" 3 , 
Co = 0.2, £| = 0.01, 6» = 0.1, p = 2.8 and D L = 4 
Gpc. Due to the relatively large values of n, £ e and £|, the 
theoretical light curve peaks at ~ 10 4 s, and it turns into 
a single line when t > 10 5 s. We see that the observed 
data points can be well fitted. 

3.5. GRB 970228 and 980326 

GRB 970228 is the first GRB that an optical counter- 
part was observed. Recently, evidence for a supernova was 
found in the reanalyzed optical and near-infrared images. 
Galama et al. (1999b) argued that the afterglow observa- 
tions are well explained by an initial power-law decay with 
a — 1.73^Q Qg, modified at later times by a type-Ic super- 
nova light curve. Here we suggest that the initial rapid 
decay of t -1 - 73 (beginning at t ~ 5.5 x 10 4 s) can be easily 
explained if a jet model was employed. This can be clearly 
seen from our calculations in Sections 3.1 — 3.4. 

GRB 980326 provides the most strong evidence for 
GRB/Supcrnova connection. Its optical afterglow has two 
distinct contributions: a power-law decaying component 
(a ~ 2.1, Groot et al. 1998) and emission from the under- 
lying supernova (Bloom et al. 1999). Similarly we suggest 
that the rapidly decaying afterglow component (beginning 
at t < 3.6 x 10 4 s) is produced by beamed GRB ejecta. 

4. Discussion and Conclusions 

Determining the existence of beaming in GRBs will be 
helpful to our understanding of the GRB "central en- 
gines" . Based on our refined dynamical model for beamed 
GRB remnants (Huang et al. 2000a), we have exam- 
ined those GRBs with rapidly decaying optical afterglows 
closely. Detailed numerical results show that afterglows 
from GRB 970228, 980326, 980519, 990123, 990510 and 
991208 can be satisfactorily fitted: the obvious break in the 
optical light curve of GRB 990123 is due to the relativistic- 
Newtonian transition of the beamed ejecta, and the rapid 
fading of afterglows from other GRBs is due to the rela- 
tively large values of £ c ,£b an d n. We thus strongly sug- 
gest that these GRBs be highly collimated. Note that in 
all cases, synchrotron radiation during the mildly relativis- 
tic and non-relativistic phases plays an important role in 
explaining the rapid decaying of optical afterglows. 

In our calculations, the £ e values distribute in a nar- 
row range, i.e., between 0.1 and 0.2. This has given some 
support to Freedman & Waxman's (1999) suggestion that 
£ has a universal value of ~ 1/3. However our results do 
not support their proposal that p has a universal value of 
2.2. Our #0 values distributed in a narrow range, ~ 0.1. 



This may provide important clues to our understanding 
of the central engine. We note that Woosley et al. (1999) 
have obtained such a small value of 9q after numerically 
studying the collapsar models of GRBs. The range of our 
n values (n ~ 1 — 1000 cm~ 3 ) indicates that some GRBs 
may be in gaseous environments, also giving some hints 
on GRB central engines. 

The jet model greatly relaxes the energy crisis for GRB 
990123 and 990510. However, we should keep in mind that 
other GRBs such as GRB 970508, 971214, 980329 and 
980703 do not have rapidly fading afterglows. They should 
not be highly collimated (Huang, Dai & Lu 2000b). Then 
the energy crisis is really a problem: GRB 971214 and 
980703 have indicated isotropic 7-ray energies of ~ 0.17 
M Q c 2 and ~ 0.06 M Q c 2 respectively! 

The rapid fading of afterglows from GRB 970228, 
980326, 980519 has also been explained as being due to the 
interaction of an isotropic blastwave with a wind environ- 
ment (Dai & Lu 1998; Chevalier & Li 1999). However the 
wind environment model could not explain the light curve 
break observed for GRB 990123 and 990510. As has been 
shown clearly in this paper, the jet model can naturally 
explain all these bursts, and should be more reasonable. 

Another possibility was proposed recently by Dai & Lu 
(1999, 2000). They suggested that the light curve break is 
due to the relativistic-Newtonian transition of an isotropic 
blastwave in a dense medium (n ~ 10 6 cm -3 ), and the 
rapidly fading afterglows can be explained as emissions 
in the non-relativistic phase. It is obviously an interesting 
proposition and should also be paid attention to. 
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